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In this (second) part of the work we present the results of numerical and qualitative analysis, 
based on a new model of the Archimedean-type interaction between dark matter and dark energy. 
The Archimedean-type force is linear in the four-gradient of the dark energy pressure and plays a 
role of self-regulator of the energy redistribution in a cosmic dark fluid. Because of the Archimedean- 
type interaction the cosmological evolution is shown to have a multistage character. Depending on 
the choice of the values of the model guiding parameters, the Universe's expansion is shown to be 
perpetually accelerated, periodic or quasiperiodic with finite number of deceleration/acceleration 
epochs. We distinguished the models, which can be definitely characterized by the inflation in the 
early Universe, by the late-time accelerated expansion and nonsingular behavior in intermediate 
epochs, and classified them with respect to a number of transition points. Transition points appear, 
when the acceleration parameter changes the sign, providing the natural partition of the Universe's 
history into epochs of accelerated and decelerated expansion. The strategy and results of numerical 
calculations are advocated by the qualitative analysis of the instantaneous phase portraits of the 
dynamic system associated with the key equation for the dark energy pressure evolution. 
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I. INTRODUCTION 

The concepts of dark energy (DE) and dark matter 
(DM) considered as two manifestations of one uni- 

fied dark fluid are basic elements of modern cos- 

mology and astrophysics. Because the dark fluid is esti- 
mated to accumulate about 95% of the Universe's energy, 
the coupling of DE with DM predestine the main features 
of the Universe evolution. In the first part of the work 
13| we introduced the so-called Archimedean-type force, 
which is linear in the four-gradient of the DE pressure 
and acts on the DM particles. This force is a relativis- 
tic generalization of the classical Archimedean force and 
belongs to the class of effective forces described in [l4| - 
[l9j . From mathematical point of view the presented 
model gives us a new self-consistent nonlinear scheme 
of interaction between two constituents of the cosmic 
medium. From physical point of view the Archimedean- 
type force is an effective redistributor of the total energy 
of the Universe between the DE and DM constituents. 
The Archimedean-type model can be also considered in 
terms of two-fluid representation of the cosmic medium 
251 ] . in which the interaction terms ±Q appear in 



the right-hand sides of separate balance equations for the 
DE and DM with opposite signs and which disappear in 
a sum, when one deals with the total balance equation. 
We formulated in [l3| the theory of interaction between 
DE and DM by using the relativistic hydrodynamics for 
dark energy and relativistic kinetics for dark matter. 

In the first paper we focused on the submodels 
with some special values of guiding parameters, which ad- 
mit exact analytical solutions of the total self-consistent 
system of master equations; in particular, we discussed 
the so-called anti-Gaussian solution. Now we consider 
the results of numerical and qualitative analysis for all 
admissible sets of guiding parameters and initial data. 



This paper is organized as follows. In Sec. II we recall 
briefly the key equations and basic formulas of the model 
with Archimedean-type coupling between dark energy 
and dark matter. In Sec. Ill the numerical results are pre- 
sented: in Sec. Ill A we explain the scheme of numerical 
analysis and details of results representation; in Sec. Ill 
B we classify the models using the number of transition 
points and discuss details of seven submodels, the per- 
petually accelerated and periodic universes being among 
them. In Sec. IV we analyze the model of Archimedean- 
type coupling in terms of dynamic system associated with 
nonlinear key equation for the DE pressure; in Sec. IV A 
we study the toy model, which relates to the autonomous 
dynamic system, find critical points and discuss two basic 
phase portraits of this system; in Sec. IV B we analyze in- 
stantaneous phase portraits of a general nonautonomous 
dynamic system in the context of numerical results pre- 
sented in Sec. Ill; in Sec. IV C we consider asymptotical 
behavior of the models. Sec.V contains discussions: in 
Sec.V A we attract an attention to the inflationary be- 
havior of the model in the early Universe; in Sec.V B we 
show that all the models under discussion manifest the 
late-time accelerated expansion; in Sec.V. C we touch on 
the coincidence problem. 



II. MASTER EQUATIONS OF THE TOY 
MODEL 



A. Dynamic parameters of the Universe 

In order to describe the Universe's evolution from the 
dynamic point of view, we find three key functions: the 
scale factor a(t), the Hubble function H(t) and the ac- 
celeration parameter —q(t) given by the standard defini- 
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tions, 



(1) 



The space-time is considered to be of the spatially 
homogeneous Friedmann-Lemaitre- Robertson- Walker 
(FLRW) type with the metric 



„3\2i 



ds 2 = dt 2 - a 2 (t)[(dx i Y + {dx 2 ) 2 + (dx 6 ) 



(2) 



When all the functions describing the state of the Uni- 
verse, S(t), depend on time through the scale factor 
<S[a(t)], it is convenient to use the new variable x and 
the following relations, associated with x: 



a(t) d , . d 

x = — 7— r, -77 = xH(x) — , t-t a = 



a(to) ' dt 



dx 



»(*o) ax 



(3) 



The last equality gives the scale factor a(t), when H(x) is 
found from the Einstein equations, which can be written 
in the form 

x—H 2 {x) = -8ttG[ P (x)+E(x)+I1(x)+P(x)} , (4) 



dx 



H 2 = 



8ttG 



(P+E) . 



(5) 



The quantities in the right-hand sides of the Einstein 
equations are interpreted as follows: p is the energy den- 
sity of the dark energy and II is its pressure; E and P are, 
respectively, the total energy density and the total pres- 
sure of the dark matter. In these terms the acceleration 
parameter —q can also be represented as 



-q(x)=l+x—logH(x)=-~ 



(6) 



B. Balance equations 

It was shown in the first part of our work that the 
separate equations of the energy balance for the DM and 
DE constituents of the dark fluid can be presented in the 
form 



x— E(x) + 3(E + P) 
dx 



x—p(x) + 3(p + U) = Q(x) 
dx 



(7) 



(8) 



The total system (DE plus DM) is considered to be con- 
served, the sum of Q and flSJ) yields the total balance 
equation 



x-r [p(x)+E(x)j + 3(p+E+U+P) = , (9) 
dx 



which is, clearly, the compatibility condition for the pair 
of equations (j4|) and ([5]). The source term Q, which was 
introduced in the first part of our work, has the form 



Q = 3x 



dx 



U(x) 



E V W P (a)(^ (10) 
(a) 



where P( a ) is a partial pressure of the DM particles of 
the sort (a) and V( a ) is the guiding parameter associated 
with the Archimedean-type force. 

The DE energy density p and DE pressure II are cou- 
pled by the linear inhomogeneous equation of state [261 ] — 

m 



p{x) = po + *U(x) + £x—IL(x) . (11) 

The state functions E(x) and P{x), which characterize 
the DM constituent, are presented by the integrals 



(a) X J ° 



} (x) e 



-A(a)V 1 +9 2 



(12) 



(13) 

where the following auxiliary quantities are introduced 



1 



F (»)( x ) = -2exp{2V (a) [n(l)-n(x)]}, 



E, 



(a) 



^(a)™(a)A( a ) 
W(a)) 



V(a) 



1(a) 



fc(B) T (a) 



p oo 

Kv(\a,)) = / dz cosh vz ■ exp [— A( a ) cosh z) 
Jo 



(14) 
(15) 

(16) 



Here i^2(A( a )) is the modified Bcssel function, k(B) is the 
Boltzmann constant, and T( a ) is t he p artial temperature 
of the particles of the sort (a) [U [33| ■ 



C. Key equation of the toy model 

Clearly, the key element of the toy model is the func- 
tion H(x); to find the DE pressure H(x) we use the key 
equation 

^x 2 U"(x)+xW{x) (4£+<r) +3(l+(j)U+3p Q =J{x), (17) 

where the prime denotes the derivative with respect to x. 
The source term J(x) — J(x, II— 11(1), II') is generally 
presented by the integral 



7W Vp [* 2 Fm(x)]' p g^Wgg 



(a) 



(18) 
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and we use them for numerical calculations. For the qual- 
itative analysis we use the explicit model expression for 
the source term 

^(x) = £*V*^^exp{V4n(l)-n(z)]}. (19) 

x^ 

Here 11(1) denotes the initial value for the DE pressure 
(let us remind that II(x=l) = U[a(t=to)]) ■ When we deal 
with massless particles, the corresponding source term 
coincides with §F5$ , if 

li = 2,, E* = E {0) , V.=V( ). (20) 

When we deal with the cold dark matter, we use the 
source term (fT9|) with 

M = 4, S. = , V, =2V (C) . (21) 

Thus, the key equation ([T7|) is the ordinary differential 
equation of the second order linear in the derivatives of 
the unknown function II" (x), II' (x) and nonlinear in the 
function II(a;) itself. 

Let us mention that, when a ^ — 1, the effective guid- 
ing parameter p$ can be introduced instead of po- In- 
deed, using a new unknown function n(x) = II(a;)— 11(1) 
one can rewrite the key equation (|17|) in the form 

£x 2 ir"(x)+xir'(x) (4£+<t) +3(l+cr)7r+3po=J r (x, tt, tt'), 

(22) 

where 

P * = po + (1 + <r)n(l) . (23) 
Clearly, the parameters /Jq and po coincide at a=— 1. 

III. NUMERICAL RESULTS 
A. About the scheme of analysis 

We analyze the Archimedean-type model as follows. 
First, we integrate the key equation (fT7|) with the source 
term (fT8|) and find the DE pressure ll(x). Then we cal- 
culate the function p(x) using (fTTj) . the functions E(x) 
and P(x) using ([12")) and (TT3]), respectively, then we find 
the Hubble function H(x) using ([5]), the acceleration pa- 
rameter —q{x) using ^ and finally the scale factor a(t) 
using ([3]). Let us make seven remarks about the scheme 
of numerical analysis. 

1. The solution for the DE pressure H(x) depends on the 
following guiding parameters: £, a, p%, V (a) , E {&) , A (a) 
and on two initial values n(l), II' (1). We fulfilled sys- 
tematic numerical calculations by varying one of these 
parameters and fixing the remaining ones, and then we 
classified the obtained plots picking out principally dif- 
ferent classes of solutions. For the numerical analysis we 
used two-component DM model, one component being 



massless (ultrarelativistic), another one being cold (non- 
rclativistic). Nevertheless, in order to illustrate the main 
details of the Universe expansion on the plots we consid- 
ered only one (leading) component of the DM. Below we 
summarize the results, emphasize the main features, but 
do not discuss details of numerical calculations. 

2. For a number of sets of the guiding parameters and 
initial data, mentioned above, the behavior of the solu- 
tion is of a singular type, i.e., the scale factor, Hubble 
function, or/and DE pressure, DE energy, etc., can take 
infinite values at a finite time moment t—t^ > t^; we do 
not discuss such solutions here. 

3. We omit the discussion of submodels for which the DM 
energy density and DM pressure grow infinitely at t — > oo 
or tend to finite (nonvanishing) asymptotic values. In 
other words, we restrict ourselves by the submodels, for 
which the DE component dominates over DM starting 
from definite time moment; this motive explains why we 
consider the parameters E^ (to ) to be unchanged in the 
presented plots. 

4. There exist a lot of sets of the guiding parameters for 
which the solutions for the DE energy and DE pressure 
are asymptotically unstable; we do not focus on them 
here. 

5. We discuss only the solutions for which the total en- 
ergy p+E is non-negative for arbitrary t, thus guarantee- 
ing that H is a real function [see (J3J)] . 

6. In order to describe the Archimedean-type model com- 
pletely we have to indicate the domain in the multidi- 
mensional space of parameters £, a, po, V( a ), Er a \, A( a ), 
n(l) and n'(l) in which the solutions are nonsingular 
and asymptotically stable; we hope to solve this problem 
in a special paper and do not discuss them here. 

7. Figures 1-7 present the results of numerical calcula- 
tions and are organized as follows. In the first panels 
labeled (a) [1 < N < 7] the plot of the DE energy p(r) is 
presented; (b) contains the plot of the DM energy density 
E(t); the DE pressure n(r) and DM pressure P{t) are 
presented in panels (c) and (d), respectively; the plot of 
the Hubble function H(t) can be found in panels (e); the 
acceleration parameter — q(r) is placed in (f), and finally, 
the plots of the scale factor a(t) are given in (g). Let us 
remind, that the parameter r = log x= log -4^ speci- 
fies the logarithmic scale in figure panels (a)-(f), and the 
value r=0 corresponds to t=to; the plots in panels (g) 
are presented in terms of logarithm of the normalized 
cosmological time t/t . 

B. Multistage character of the Universe's evolution 

The analysis of the pictures presented in 
Figs.l(f),2(f),3(f),4(f),5(f),6(f), and 7(f) shows ex- 
plicitly that for a wide choice of sets of the parameters 

v, Po, V( a ), £( a ), A( a ), n(l) and n'(l) the function 
—q[t) is not monotonic and can change the sign. We 
indicate the points, in which the function —q(r) changes 
the sign, as transition points, which divide the history 
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of the Universe into deceleration/acceleration epochs; 

(k) 

we denote these points as 7V trans v The epochs with 
accelerated expansion can give way to the epochs 
with deceleration, the number of changes being finite 
or infinite. When the (continuous) curve —q(r) has 
extrema, it is natural to divide the corresponding epochs 
into eras, the number of eras relates to the number of 
maximums and minimums. Thus, depending on the 
number of transition points (related to the number of 
(simple) roots of the equation q(r) = 0), we divide below 
all the models into three classes and into a number of 
sub-classes. 



1. Perpetually accelerated universe (no transition points) 

The first class of models relates to the case when the 
function — q(r) is non-negative for r > 0. In [l3T ] we 
discussed the so-called anti-Gaussian model: exact so- 
lutions in this model were presented in terms of ele- 
mentary functions at a——l. In particular, the accel- 

3V 

eration parameter obtained there, —q(t)=l+ WlT Q(fl t tyi , 
is a monotonic positively defined function exceeding the 
value — (?(dS) = lj the acceleration parameter for the de 
Sitter models. In other words, there are no transition 
points, there is no splitting of the history into epochs 
and eras, when we deal with an anti-Gaussian universe. 
Another example presented in [l3j describes superexpo- 
ncntial expansion with the acceleration parameter given 

by -q (t )=1+,/ «fiL If- Z=. — °Jl. The plot of 

* W V W) cosh 2 |yi27rGpo(*-to)J 

this function starts with —q(ta)=l, reaches the maxi- 
mum -(? (max )=l+y / ^fy at t=t o + l °^l^\ and tends 

asymptotically to — g(oo)=l. In the context of our ter- 
minology, we also deal with perpetually accelerated uni- 
verse, but its history is now divided into two eras. Here 
we discuss an alternative example of a perpetually accel- 
erated universe with a ^ — 1, studied numerically. Typi- 
cal pictures illustrating this class of models are given in 
Fig.l. 

The solid curve in Fig. 1(f) touches the line q = in 
one point, other curves lie above q — 0, thus there are no 
transition points. In this class of models the acceleration 
parameter —q as a function of time is non-negative for 
arbitrary r: this means that the accelerated expansion of 
the Universe is perpetual. For the model under discussion 
the DE energy density p is (perpetually) non-negative 
[see Fig. 1(a)], the DE pressure II is nonpositive. The 
acceleration parameter, the DE energy density and pres- 
sure, and the Hubble function change nonmonotonically 
with time, thus, the statcfindcr parameter r = ^t^t, 
introduced in 34], is necessary to describe the details 
of the Universe's evolution. We do not focus now on 
these fine details of the model and their physical inter- 
pretation, nevertheless, we have to stress, that extrema 
of the function —q(r) divide naturally the history of the 
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FIG. 1. Perpetually accelerated universe: the acceleration 
parameter —q(r) is non-negative for arbitrary r. Since there 
are no transition points, in which — q(r) changes the sign, 
the history of this universe includes only one epoch; since 
at least five extrema can be recognized in Fig. 1(f), the his- 
tory consists of six or more eras. Here the dark solid line re- 
lates to the following set of the parameters: £=0.35, a—— 0.99 
(i.e., ~ 0.06 > 0), V (0 )=l, £ (0) =0.0205, A (0) =l (i.e., 

the DM is initially medium-relativistic), p*=0.333 • 10 -4 , and 

n'(i)=-i. 



Universe into the corresponding eras, thus describing the 
multistage character of its evolution. 

The history of a perpetually accelerated universe in- 
cludes only one epoch, but consists at least of six eras. 
Indeed, one can clearly identify at least three maxi- 
mums and two minimums of the curve —q(r) which di- 
vide the acceleration epoch into six eras. The first era 
(t < 10°) can be indicated as the era of a first superac- 
celeration. During this era the DE energy - density, p(r), 
grows monotonically [see Fig. 1(a)], the DM energy den- 
sity, E(t), decreases [see Fig. 1(b)], but at the beginning 
of the era DM dominates over the DE from the point of 
view of energy contribution. The second era finishes with 
slumps of the DE energy density, of the Hubble function 
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and of the acceleration parameter. The DE pressure, 
which reached its global maximum at this moment, be- 
comes, clearly, a locomotive power of further universe 
evolution. The third era can also be indicated as the 
era of super-acceleration, which is characterized by the 
growth of DE energy density and the Hubble parame- 
ter. Fourth, fifth, etc. eras can be considered in terms 
of relaxation to the state with asymptotically constant 
positive values — g^, i?oo, Poo and negative IIoo. The 
behavior of the DM state functions E(t) and P{t) dur- 
ing the second, third, etc. eras can not be recognized in 
the plots in Fig. 1(b) and Fig. 1(d), that is why we visual- 
ize their behavior specially in the discussion below using 
more appropriate scale (see Fig. 11 and Fig. 12). 

2. Periodic models (infinite number of transition points) 

We indicate the model by the term periodic when the 
equation q(r) = has an infinite number of roots and the 
history of the Universe splits into infinite number of iden- 
tical epochs with accelerated and decelerated expansion. 
Clearly, it is the case opposite to the first one (perpetu- 
ally accelerated universe) . Figure 2 illustrates this termi- 
nology: the DE energy density and DE pressure, as well 
as the Hubble function and acceleration parameter oscil- 
late with fixed frequency and amplitude (we should keep 
in mind that the time scale in Fig. 2 (a)-(f) is logarithmic, 
so that the frequency seems to be visually increasing) . In 
the model presented in Fig. 2 the DE energy density p is 
non-negative, the averaged over time (mean) value of the 
DE pressure < n > is negative, and the sum p + II is 
positive. The plots in the Fig. 2(e) and Fig.(2)(f) demon- 
strate that the stable oscillatory regime starts in the Uni- 
verse after the second transition point. The first epoch is 
characterized by the decelerated expansion and consists 
of three eras; at the end of the first era (r ~ 10°) the 
DM contribution into the total energy density becomes 
negligible in comparison with the DE contribution. The 
second epoch relates to the accelerated expansion and can 
be split into two eras, the maximal value of the accelera- 
tion parameter being about 3 times smaller than for the 
subsequent epochs. The scale factor a{t) oscillates near 
the de Sitter curves [see Fig. 2(g)]. The periodic model 
can be obtained as a limiting case, when the parameter 
3£+<7 tends to its critical value 3£+cr — > 0+o- 



3. Universe evolution with a finite number of transition 
points 

Let us focus on the models for which the equation 
q(r) = has a finite number of simple roots, m > 1. 
This means that during the Universe evolution the epochs 
with accelerated expansion are changed m times by the 
deceleration epochs, and vice versa. When the roots are 
not simple, special attention to such models is necessary, 
since the existence of the root does not guarantee the 
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FIG. 2. Periodic universe with infinite number of acceler- 
ation/deceleration epochs. This model is a limiting case, 
when the parameter 3£+er tends to its critical value 0+o ■ For 
this model not only the sum p+E, but the DE energy den- 
sity p itself remain non-negative. Oscillations in the plots of 
the Hubble function H(t) and of the scale factor a(t) are 
the most explicit (see dark solid line), when the parame- 
ters of the model are the following: f=0.1, cr=-0. 299999 
(i.e., =i 10" 6 > 0), V( )=l, £ ( o)=0.0205, A (0) =l, 

p„=0.333 • 1CT 4 , and n'(l)=0.01. 

change of the sign of the function —q{r) (e.g., when the 
root is double, the curve —q(r) touches the line q = 0, 
and this root does not give a transition point). The 
class of models with m transition points can be natu- 
rally divided into subclasses with m=l, m=2, etc.; every 
subclass can be divided into sub-subclasses according to 
physical and geometrical motivation. Let us focus on the 
basic ones. 

(i) One transition point. 

Figure 3 illustrates the subclass of models, for which 
the plot of the acceleration parameter —q(r) (see 
Fig. 3(f)) can be visualized as a deformed Heaviside step- 
function. The change of the deceleration epoch by the ac- 
celeration epoch takes place only once in a narrow time 
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FIG. 3. The model with one transition point and —q(r) curve 
of the Heaviside step-function type. The first (deceleration) 
epoch and the second (acceleration) epoch are not divided 
into eras. The epoch with the accelerated expansion is of the 
de Sitter-type, characterized by p+IT=0. The sum p(t)+E(t) 
is non-negative for arbitrary r > 0. The parameters of the 
model are the following: £=0.1, a— 50, V(o)=l, -E( )=0.0205, 
A (0) =l, p„=0.333 • 10~ 4 , and lT(l)=-5. 



period, and the plot looks like a typical picture for a 
phase transition. Taking into account the plot of the 
Hubble function [see Fig. 3(e)] one can conclude that the 
Universe's evolution starts and finishes on quasi de Sitter 
asymptotes with different Hubble constants H ^ H^. 
The functions p(r), n(r), related to the DE component 
of the dark fluid, monotonically increase, and the func- 
tions E(t), P(t), related to the DM, monotonically de- 
crease, the sum p(t)+E(t) being non-negative for arbi- 
trary t > 0. The second epoch, characterized by the 
accelerated expansion, looks like the de Sitter-type stage 
with p + n = and vanishing E and P. The first and 
second epochs are not divided into eras in the framework 
of this model. 

Figure 4 illustrates an alternative subclass of models 
with one transition point. In contrast to the first ex- 



FIG. 4. The model with one transition point and two eras 
in both decelerated and accelerated epochs. The model de- 
scribes the transition from one de Sitter-type universe to an- 
other one with different Hubble constants Ho =fc Hoc- In 
contrast to DM energy density, the DM pressure grows in 
the early universe, reaches the maximum and then decreases 
quickly in the deceleration epoch. The parameters of the 
model are the following: £=0.1, cr=0, V( )=l, S(o)=0.0205, 
A (0) =l, p„=0.333 • 10~ 4 , and n'(l)=-17. 



ample, the functions p(r), n(r), H(t), P(t) and — q(r) 
have extrema, so that one can pick out two eras in both 
epochs of the Universe history. This model also describes 
the transition from one de Sitter-type universe to an- 
other one with different Hubble constants Ho ^ i?oo- It 
is important to stress that the DM pressure P(t) is not 
monotonic now: in contrast to the DM energy density, 
the DM pressure grows in the early Universe, reaches the 
maximum and then decreases quickly in the deceleration 
epoch. 

(ii) Two transition points. 

The example of the model with two transition points is 
given in Fig. 5. The behavior of the functions p(r), n(r), 
E(t) 7 and H{t) is similar to the one, given in Fig.l, but 
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FIG. 5. The model with two transition points, two accelera- 
tion epochs and one deceleration epoch. The plots for the DE 
energy density, DE pressure, and Hubble functions are analo- 
gous to the curves given in Fig.l. The DM pressure has a max- 
imum at the end of the first era of the first acceleration epoch 
in analogy with the curve given by Fig. 4(d). The parame- 
ters of the model are the following: £=0.1, a=— 0.08, V(o)=l, 
£ (0) =0.0205, A (0) =l, p«=0.333 ■ 10^ 4 , and n'(l)=-15. 



there are new details in the panels (d) and (f). First, the 
DM pressure reaches maximum at the end of the first era 
of the first acceleration epoch. The second new detail 
is that the deceleration epoch appears, which separates 
two acceleration epochs with different asymptotic values 
of the Hubble function. 

(Hi) Three transition points. 

The model which displays two epochs of decelerated 
expansion and two epochs of accelerated evolution, is pre- 
sented by Fig. 6. The graph of the function —q(r) is of 
the so-called AT type. The start of the Universe's ex- 
pansion relates to the deceleration epoch, which is then 
replaced by a short acceleration epoch as in the case of 
the periodic model [compare Fig. 6(f) with Fig. 2(f)]. The 
Universe's evolution during the second epoch of decel- 



FIG. 6. The model with two epochs of decelerated expansion 
and two epochs of accelerated evolution. The plot of the func- 
tion — q(r) is of the so-called Af type. The second accelerated 
epoch is of the de Sitter type. The parameters of the model 
are the following: £=0.1, a=l, V (0 )=l, -E (0 )=0.0205, A (0) =l, 
p,=0.333 ■ 10" 4 , and n'(l)=-5. 



erated expansion and the subsequent second accelerated 
epoch is analogous to the model with one transition point 
(compare Fig. 6(f) with Fig. 3(f)). Clearly, the second (fi- 
nal) accelerated epoch is of the de Sitter type. 

(iv) Quasiperiodic evolutionary models. 

Figure 7 illustrates the behavior of the model with a 
big but finite number of transition points. The first, 
second, and third epochs of deceleration, as well as the 
first, second, and third epochs of accelerated expansion 
in this model have the same character as for the peri- 
odic model [compare with Fig. 2(f)]. Nevertheless, start- 
ing from some transition point with the number (k), i.e., 

when t > t^rans)' ^ ne curve ~ a ( T ) remains above the line 
q=0, thus, later the Universe's expansion is accelerated. 
Clearly this model is intermediate between the model 
with two transition points [see Fig. 5] and periodic model 
[see Fig. 2]. The behavior of the functions p(r), n(r), 
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FIG. 7. The intermediate model with big but finite number 
of transition points. The oscillatory regime is produced by 
the DE component of the dark fluid, and it is switched on 
at the moment when the DM energy density and DM pres- 
sure become vanishing. Late-time expansion is, evidently, ac- 
celerated. The dark solid lines relate to the following pa- 



rameters: £=0.1, 



£ {0) =0.0205, A 



(o) 



ct=-0.29 (i.e.. 
=1, p»=0.333 ■ 



3C+cr ~ 0.01 > 
10 -4 , and n'(l): 



0), V, 
=-1. 



(0) = 



H(t) and —q(r) is quasiperiodic: the amplitudes of their 
oscillations decrease asymptotically. Late-time behavior 
in average is of the de Sitter type with positive DE en- 
ergy density and negative DE pressure. The interest in 
such models has been renewed after publication of the 
paper [35| . 



IV. QUALITATIVE ANALYSIS 

In Sec. Ill we classified the solutions obtained nu- 
merically according to the number of transition points, 
in which the decelerated expansion of the Universe is 
changed by the accelerated expansion and vice versa. In 
order to explain the results qualitatively, below we ana- 
lyze the key equations by two methods: constructing the 



phase portraits and using asymptotic analysis. First of 
all, we consider the example of the model with a = — 1 , 
for which the key equation can be reduced to the au- 
tonomous dynamic system; we find critical points and 
draw the corresponding phase portraits. On the basis 
of these results we study nonautonomous dynamic sys- 
tems with (7^-1 taking into account that critical points 
of this nonautonomous system remain of the same type 
as for the case o= — 1, but they drift with time on the 
phase plane. Finally, we analyze three asymptotic limits: 
LT — » oo, II — > and II — > const ^ 0, which allow us to 
simplify the key equation significantly. 



A. Dynamic system associated with the key 
equation 

Let us introduce new variables X(r), Y(t) and r as 
follows 

X = 3^% x~'"- 1 exp{-V*[0(z) - 0(1)]} , 



Y 



V* 



dn 



(/i+l) dx 



— , t= (/i+l) log x. (24) 



Then the key equation of the second order ([TT]) with the 
source term (|19p can be rewritten as a pair of dynamic 
equations of the first order: 



dX 
~d7 



-X(l + Y) 



r+ lo g |%±Ux 



where 



.4 



£(M + 1) 2 



B 



(25) 
(26) 



The quantity X(t) is non- negative due to its definition 
(PMl) . X tends to zero in two different cases: first, when 
n — > oo, or second, asymptotically at x — > cxj, if n(cx)) is 
finite. The initial value for the parameter r, t—0 relates 
to the moment t—tg and, respectively, x=l, then initial 
data for the integral curves in terms of these variables 
are the following: 



X(0) = 



E^ V* 



Y(0) 



V, 



(/i+l) 



O'(l). 



(27) 



Generally, this dynamic system is nonautonomous be- 
cause of the term linear in r in the second equation. 
Nevertheless, there is a case er=— 1 when the system |25|) 
becomes autonomous. Let us consider this case in more 
detail. 
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Autonomous dynamic system at a = — 1 



The autonomous dynamic system 



dX 



= -X [1 + Y] 



dY 



= -A + Y[X-B] (28) 



contains only two control parameters: A and B, the pa- 
rameter A being now positive. The new dimensionless 
parameter 



B = (3g-l)( M +l) 
A 3 Po V* 



(29) 



plays an interesting role in further analysis. Physically 
appropriate integral curves should be situated in the do- 
main X > on the plane XOY. We consider the integral 
curves, which penetrate into the domain X < 0, as non- 
admissible and omit the models with the corresponding 
set of guiding parameters. The qualitative analysis of the 
two-parameter autonomous dynamic system (|2"51) demon- 
strates the following interesting features. 



2. Critical points 

When A ^ B and B ^ 0, there are two critical points 
associated with this dynamic system. The first critical 
point 



*(D=0, Y (1) = -- = 



(i) 



.4 
B 



(30) 



is the saddle one, when 9 < 1, and is the stable nodal 
point, when 9 > 1. The second critical point is 



X {2) =B-A = 



-1) 



3poV* 



Y m = -i 



(31) 



and it coincides with the first one, when B=A, i.e., when 
9=1. Generally, the second critical point can be 
(i) the saddle point, when 9 > 1 (the first critical point 
is the nodal one); 

> 1-9 > 



3poV, 



(ii) the stable nodal point, when ^^ +1 yi 
(the first critical point is the saddle one); 



3p V» 



(Hi) the stable focus, when < 4 ^^ +1 ^ 
critical point is the saddle one); or 
(iv) the degenerated node, when < j^^fy—l 
first critical point is the saddle one). 



< 1-6 (the first 
(the 



3. Two phase portraits of the autonomous dynamic system 

A typical phase portrait of the autonomous dynamic 
system with a stable node and a saddle point is pre- 
sented in Fig. 8. Integral curves in the northern and 
western domains, which are cut out by the separatrices of 
the saddle point with coordinates X( 2 )=0.2125, Y( 2 )=—l, 
tend asymptotically to the nodal point with Xm=0, 




FIG. 8. Phase portrait of the autonomous dynamic system 
with stable node at X(!)=0, Yn\=— 0.15 and saddle point at 



X {2) =0.2125, 



Y(2)=— 1. Integral curves in the northern and 
western domains, which are cut out by the separatrices of the 
saddle point, tend asymptotically to the node. Integral curves 
in the southern and eastern domains of this saddle point can 
be asymptotically characterized by infinite values of X and/or 
Y. 



Ym=— 0.15. All the integral curves in the southern and 
eastern domains of this saddle point can be asymptoti- 
cally characterized by infinite values of X and/or Y . 

A typical phase portrait of the autonomous dynamic 
system with a focus and a saddle point is given in Fig. 9. 
The domain of quasiperiodic motion is situated in the 
south-western domain, which is cut out by the separa- 
trices of the saddle point at Xm=0, Fm=0.25. This 
quasiperiodic motion is associated with the presence of 
the focus at X@)= — 0.59375, Y(2)=— 1. Curves walking 
in three other domains are asymptotically unstable. 



B. On the behavior of the nonautonomous 
dynamic system with a > — 1 

When cr=— 1 and the dynamic system is autonomous, 
there are two typical phase portraits, given by Fig. 8 and 
Fig. 9. The qualitative difference between them is pre- 
determined by specific combinations of the guiding pa- 
rameters £, po, and V*: any changes on the phase space 
are associated with variations of these guiding parame- 
ters. Typical phase portraits with fixed guiding parame- 
ters remain unchanged with time, and the integral curves 
can be prolonged till r — > oo. When a > — 1, the dy- 
namic system is nonautonomous, nevertheless, one can 
describe such a model qualitatively in terms of instanta- 
neous phase portrait distortion, in terms of (quasi) critical 
points drift, and their transformations. One can stress 
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f (X) = X - p log X 
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FIG. 9. Phase portrait of the autonomous dynamic system 
with stable focus and saddle point. The domain of quasiperi- 
odic motion is situated in the south western domain, which is 
cut out by the separatrices of the saddle point at Xnj=0, 



yn)=0.1875. The focus is situated at 



Y, 



X (2 ) = 



0.59375, 



(2) = 



-1. 



that in this case a restructuring of the instantaneous 
phase portraits can take place because of time growth, 
when the guiding parameters remain unchanged. 



FIG. 10. The plot illustrates the appearance of a pair 
of (quasi)critical points, when the cosmological time in- 
creases. The function f(X) = X—f3logX has the minimum 
/(min)— ft log j at X=P; the point of minimum divides the plot 
of this function into two branches: first of them is situated 
at < X < ft, the second one relates to the interval X > ft. 
Clearly, Eq. (|32[) has no solutions, when K,(tq) < /31og4; it 
has only one (double) root, when /C(ro)=/31og -|; there are 
two different roots, when K.(tq) > /Slog §, the corresponding 
points indicated as Mr*) and Mr»*y The curve on the plot 
corresponds to the value /3 = 1. 



1. Critical points drift and transformation 

Let us consider a close vicinity of the point t=to and 
put the value To to the right-hand side of the second 
equation (1251) . In the small vicinity of To (t=to+t;, and 
?/ << To is a small deviation) we can consider (quasi) 
critical points in analogy with true critical points, stud- 
ied in the previous case. The first (quasi) critical point 
relates to F( 1 )=oo, when 3£+<7 < 0, and 

Yn)=— oo, when 3£+cr > 0. In both cases these points 
disappear from the admissible part of the dynamic plane 
X0Y. Other (quasi) critical points relate to the value 
Y(2)=— 1, and the corresponding values X can be found 
from the transcendent equation 



X-p\ogX = K{T ) 



3(l+o-) 
£(M+1) 2 



>0, 



(32) 



K(t )=/3 



T + log 



gQi+1) (3g+<7)(/x+l)-3p5V. 



3(1+0-) 

(33) 

The plot of the function f(X) = X— filogX is pre- 
sented in Fig. 10. The function f(X) has the minimum 
/(mm) =/8 l°g § at X—f3, which divides the corresponding 
curve into two branches: the first of them is situated at 



< X < f3 (left branch), and the second one relates to 
the interval X > f3 (right branch). Clearly, Eq. (|32|) has 
no solutions, when K{tq) < (3 log % or equivalently when 
r o < T (crit) i where the critical time moment is 



' (crit) : 



d | plV* ( M +l)(3e+a) [ (»+l)V,E, 



1+CT 



3(l+a) 



3(l+<r) 



(34) 



Eq. (|3^|) has only one (double) root, when /C(To)=/31og 4, 
or T =T( crit ); there are two different roots, when JC(t ) > 
/31og % or equivalently To > T( crit ). In other words, there 
are no (quasi) critical points at To < T( crit ), one point ap- 
pears at To=T( cr it), and two (quasi)critical points appear 
in the instantaneous phase portrait, when To > T( cr j t ). 

The first (quasi) critical point indicated as M(*)(to) 
drifts to infinity, since the value for which the right 
branch of the curve happens to be crossed by the hori- 
zontal straight line, increases with To. In order to recog- 
nize the type of this (quasi) critical point, let us find the 
roots of the characteristic equation associated with the 
dynamic system in the vicinity of the point with coordi- 
nates (X(*), -1), 



Ai,2 - - 



X^—B ± 



^/(X ( , ) -i?) 2 +4(X (t) -/3)J . (35) 



At the critical moment of time the equality X^=f3 takes 
place, and we deal with the degenerated node; when 
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> j3, there are two real roots of different signs, thus 
the (quasi) critical point is the saddle one. 

The second (quasi)critical point Mu*\ drifts to the 
point (0,-1) when To increases, since the value Xu^ < 
j3, for which the left branch of the curve crosses the hori- 
zontal straight line, tends to zero. The type of the second 
(quasi)critical point depends on the value To and on the 
relation between parameters B and j3. 
(i) When B > 1+/3, or cquivalently 



a > 



/'- 



r+£(M+l) 



(36) 



the discriminant in (|35p is always positive, the roots of 
characteristic equations are real and negative, thus the 
second (quasi) critical point is a stable node for each To > 
T (crit)- We deal with the phase portrait, which contains 
a stable node Mu*\ and a saddle point Muy 
(ii) When B < 1+/3 we face the reconstruction of the 
phase portraits with time. Indeed, one can fix two mo- 



ments of time, Tn and r!f J , defined as 



(2) 



XuAt^) = B-2 + 2y/l + P-B : 



X { , h (T^) = B — 2 — 2^1 + p-B. 



(37) 



for which the discriminant in (1351) takes zero value. Since 



X 



(**) 



tends to zero with time, one obtains three subse- 



quent situations: first, when T( crit ) < To < Tq the roots 
are real, i.e., M(*„) is a node; second, when Tq < to < 

(2) 

Tg , the roots are complex numbers, i.e., Mu^ is a fo- 
cus; third, when T( 2 ) < to the roots are real again, i.e., 
M(**) is a node. One can guarantee, that the nodes and 
focus are stable, if (3 < B < 1 + /?. 

(iii)lf B < /?, evidently, there exists a moment Tq , when 
the decreasing quantity X^^(tq) < f3 becomes equal to 
the parameter B. At this moment A1.2 = H\//3—B, i.e., 
the focus transforms into a center. It is possible, when 



1 < a < 3 



(M-2) 



(38) 



^ and Tn* 1 the discriminant 



Finally, at the moments Tq 1; and Tq 2 ^ 
in (|35p vanishes, i.e., the corresponding points should be 
classified as degenerated nodes. 



2. Phase portrait distortion and transformation 

Let us consider the plots in 
Fig.l(c),2(c),3(c),4(c),5(c),6(c), and 7(c), obtained 
numerically, from the point of view of qualitative anal- 
ysis of the models with a > — 1. In fact we take two 
typical phase portraits: for a node and a saddle point 
[Fig. 8], and for a focus and saddle point [Fig. 9], and 
provide the changes of two types. First, we move the 



(quasi)critical points horizontally, thus distorting the 
instantaneous phase portrait; second, we replace the 
node by focus and vice-versa, thus modeling transitions 
in the instantaneous phase portraits with time. 

1. Figurel(c): perpetually accelerated universe. 

For the given values of the guiding parameters we have 
a < 3+4£, i.e., the inequality B < 1+/3 is valid [see the 
discussion of Eq. (|3"T|) above]. The zone t < 10° can be 
indicated as precritical, the function II(t), the DE pres- 
sure, decreases there, since the initial value LT'(1)= — 1 is 
negative. Then, as it was qualitatively predicted, the fo- 
cus appears on the instantaneous phase portrait, so that 
the DE pressure enters the zone of damped quasioscilla- 
tions. 

2. Figure 2(c): periodic universe. 

For this model, clearly, both inequalities given by (|3"8")l 
are valid; in addition B < 1+/3, thus it is a typical model 
with center [see the discussion od Eq. (|3"5|) above]. The 
presence of a center on the phase portrait explains the 
periodic graph of the DE pressure II(t) (the graph is 
presented in a logarithmic scale). 

3. Figure 3(c): first model with one transition point. 
There, clearly, the inequality (|3l)]) is valid, i.e., B > 1+/3 
and we deal with the absence of focuses. Instantaneous 
phase portraits contain stable nodes only, the graph for 
the DE pressure has no quasioscillations. 

4. Figure 4(c): second model with one transition point. 

5. Figure 5(c): model with two transition points. 
There, in both cases, we deal with the situation similar 
to the one in Fig. 1(c). 

6. Figure 6(c): model with three transition points. 

The situation is analogous to the case given by Fig. 3(c). 

7. Figure 7(c): quasiperiodic universe. Again we deal 
with the set of guiding parameters, for which B < 1+/3, 
and thus the stable focus appears in the instantaneous 
phase portrait of the dynamic system; quasioscillations 
of the DE pressure are typical for this model. 



C. Asymptotic behavior of the integral curves 

Speaking about nonrelativistic and ultrarelativistic 
models, we distinguish between initially and effectively 
hot or cold states of the DM particles. When the param- 

are much bigger or much smaller 



eters A( a ) = 



,T(a 



than 1, we deal with initially nonrelativistic or ultrarel- 
ativistic DM, respectively. In the course of the Universe 
expansion an initially cold DM can become effectively hot 
(effectively ultrarelativistic), when the term (x) — > 00 
[sec (|14p]. Such terminology is advocated by the fact, 
that at F^(x) —> 00 one obtains < q 2 > Fu\(x) >> 1, 
and the averaged kinetic energy of a DM particle be- 
comes much bigger than the corresponding rest energy, 
even if initially at t=to (F(l)=l) the inequality had the 
form < q 2 ><< 1. Analogously, when Fm(x) — > 0, 
an initially hot DM converts into the effectively nonrel- 
ativistic, because the rest energy dominates over the ki- 
netic energy of a DM particle. A new interesting case 
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can appear due to the Archimedean-type force, when 
F(a)( x ) ~~ ^ const ^ 0, despite the fact that a(t) — > oo. 
In this case one deals with the asymptotic fixation of the 
degree of relativism of the DM particles instead of stan- 
dard effective cooling. Let us consider these asymptotic 
cases and analyze the model behavior for arbitrary values 
of the parameters A( a ). 

1. Asymptotic regime F( a )(x) — > oo 
In the absence of the Archimedean- type force, i.e., 



at V, 



(a) 



0, the function Fu\(x) 



tends to zero 



asymptotically. When V( a ) ^ 0, the function F^(x) can 
tend to infinity if the DE pressure takes negative infinite 
value, n(x) — > — oo. In this asymptotic regime, when 



the integrals in ([11]), (fT3]), (fTH]) can be 



easily calculated yielding 



J(x) 



U'(x) 



X^(a) V (a) J (a) e 

(a) 



v (a) [n(i)-n»] 



E(x)=3P(x) = —Y,E {a) I (ll) e 



v (a) [n(i)-n(z)] 



(a) 



where the parameter J( a ) is given by 



r (a) = 2A^[Af a) +3A (a) +3]e- 



(39) 



(40) 



(41) 



Let us indicate the maximal parameter from the set 
{V( a )} as v, the corresponding parameter as E and 
/(a) as /. Then one obtains that, when a > — 1 and 
x — » oo, the leading order term in the left-hand side of 
the key equation (IT7|) is 3(l+er)II and the reduced equa- 
tion for II(x) is of the form 



U'(x) = 



3(1 



vEI 



J n(a;)exp{^n(a;)} . 



(42) 



This means that at II — s- — oo the derivatives II', LI", 
etc., asymptotically vanish, and the behavior of the DE 
pressure II(a;) is described by the function inverse to the 
integral exponent i?i _1 (II), where Ei(U) = J ^-e~ vY1 . 
The decomposition of this function at II — > — oo gives 
the transcendent equation 



exp {— vll(x)} 



3(1 + a) 
4EI ' 



(43) 



When a = — 1, the asymptotic law is H{x) — » —-logx in 
accordance with the exact solution of the anti-Gaussian 
type, studied in the first part of the work [lijj ]. The 
asymptotic behavior of the DE energy density is given 
by the formula p{x) ~ oTI [see (fTTj) ] and depends essen- 
tially on the sign of the parameter a: it remains positive 
at II(x) — > — oo, when — 1 < a < 0. Concerning the DM 
energy density E(x) and DE pressure P{x), these quan- 
tities grow exponentially; we do not consider the Hubble 
function H(x) and the scale factor a(t) for such models, 
since there is no observational data, which could confirm 
such behavior. 



2. Asymptotic regime F^(x) — !> 

This asymptotic regime can be realized, when the DE 
pressure takes arbitrary positive or finite negative value 
at x — > oo. The asymptotic behavior of this type relates 
to the vanishing source term 

311' 

J(x) ->■ — AAexp{2V (a) [n(l)-n(.x)]} , 



(a) 



(a)fc(B)^(a)V( a) 



^3(A (a) ) 
*2(A (a) ) 



and vanishing DM energy density 

(a) 



(44) 



(45) 



Thus, the DE pressure satisfies now the Euler equation 
£x 2 ll"{x) + xll'(x) (4£ + < 7)+3(l + cr)n + 3 / o = 0, (46) 

which is studied in detail in Sec. Ill A of the first part 
of our work [jjjj . In particular, in Sec. Ill A 3, we have 
shown that there exist quasiperiodic oscillations near the 
asymptote II(x) = — which are manifestly presented 
in Fig. 1(c), Fig. 4(c), Fig.5(c), Fig.7(c). The aperiodic so- 
lutions studied in Sec. Ill A 1, and Sec. Ill A 2 [l3[ clearly 
appear in Figs. 3(c) and 6(c). 

3. Asymptotic regime F^(x) — > const ^ 

Let us suppose for simplicity that V( a )=i^ for all types 
of the DM particles, and let the corresponding function 
F(x) tend to a nonvanishing constant value Fqo, when 
x — ^ oo. This means that asymptotic behavior of the DE 
pressure is the following 

U(x) ->■ n(l) - - log x - — log F^ . (47) 
v 2v 

Then the initial data require that i 7 ^ = 1, and the source 
term (fTB|) takes asymptotically the following form 



J(x) 



(48) 



(a) 



i.e., it decreases as x~ 3 at x — > oo. The left-hand side of 
Eq. ([T7]) vanishes after the substitution of (|4"T|) . if a = —1 
and 3pois = 3£— 1, i.e., the asymptotic regime exists at 
special conditions only. The DE and DM energy densities 
behave in this case as 

p(x) -> ^logx, 



E(x) 



(a)fc(B)^ (a) 



(a) 



Ha) 



^3(A (a) ) 
^2(A (a) )' 



(49) 
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The Hubble function, scale factor and acceleration pa- 
rameters are given by the following formulas 



/ 87rG f 2ttG 2 i 

H(x) ->• \j-^-\ogx, a(t) -> a(i ) cxp j i f> , 



9(0 -► 1 + 



3;y 



47rGt 2 ' 



(50) 



which, clearly, describe the anti-Gaussian solution stud- 
ied in the first part of the work (13j . 



a(t)/a(0) 



E(t), P(x) 




V. DISCUSSION 

The Archimedean-type model can describe qualita- 
tively two cornerstones of the Universe's evolution: the 
inflation in the early Universe and the late-time acceler- 
ated expansion. Also, this model offers a variant expla- 
nation of the so-called coincidence problem, and it gives 
some motivation for the specific division of the Universe's 
history into epochs and eras. Let us consider these ques- 
tions in more detail. 



A. Inflation in the early Universe 

When r < 1 Fig.l(a), 2(a),3(a),4(a),5(a),6(a), and 7(a) 
display the first inflationary-type epoch: one can reveal 
a quick growth of the DE energy density p(r). During 
this period of the Universe evolution the scale factor a(i) 
increases essentially: in order to visualize this inflation- 
ary effect we presented additionally the fragment of the 
plot a(t) in the natural scale t instead of the logarithmic 
one [see Fig. 11 (a)]. The steepness of the plot of a(t) is 
predetermined by the combination of the guiding param- 
eters; there are a lot of possibilities to fit the observed 
curve by the curve predicted in the framework of the 
Archimedean-type model. In the context of expansion 
the growth of the DE energy density p(r) in the early 
Universe relates to the heating of the DE component of 
the dark fluid. Clearly [compare, e.g., the panels (a) 
and (b) in Figs. 1-7] the sum p+E remains non- negative 
for the whole interval of time. This feature provides the 
H 2 to be non-negative, which guarantees that H is a 
real function. At the same time the DM energy density 
E(t) [see the panels (b)] decreases monotonically; never- 
theless, the rate of its effective cooling differs essentially 
from the rate given by the standard powerlaw function 
(1/a 4 for the ultrarelativistic DM or 1/a 3 for the cold 
DM). Principally new detail can be found in Fig. 4(d) and 
Fig. 5(d): the DM pressure is described by a nonmono- 
tonic function; it grows quickly, reaches a maximum, and 
then decreases rapidly. 



FIG. 11. In the left panel a fragment of the plot is presented, 
which illustrates the inflationary-type growth of the scale fac- 
tor a(t) for the early Universe, given in terms of cosmological 
time t. In the right panel one can see a typical late-time be- 
havior of E(t) and -P(t); this rescaled fragment illustrates 
the features of the DM evolution at r > 1, which cannot be 
recognized on the plots given by Fig.l(b,d)-Fig.7(b,d). 



B. Late-time accelerated expansion 

When r — > oo, the curves —q(r), which present 
the evolution of the acceleration parameter in 
Fig.l(f),2(f),3(f),4(f),5(f),6(f), and 7(f) tend to the 
horizontal asymptote, the value — g(oo) being positive. 
Taking into account the behavior of the plots of the 
Hubble functions displayed in Fig. 1(e), Fig. 3(e)-Fig. 7(e), 
one can conclude that H (r — > oo) — > -f/oo, he., the Hub- 
ble function also tends to the positive constant value, 
when t — y oo. This means, first, that a typical behavior 
of the Universe with Archimedean-type interaction 
between DE and DM is characterized by the late-time 
accelerated expansion; second, that asymptotically the 
Archimedean-type model converts into the quasi-de 
Sitter one. 

The DM state functions E(t) and P(t) decrease 
rapidly, and for r > 1 we should change the scale in 
order to visualize the behavior of these functions. A typ- 
ical late-time behavior of E(t) and P(t) is presented in 
the right panel of Fig.ll. In order to describe the effec- 
tive equation of state of dark matter in the framework of 
the Archimedean-type model, we calculated numerically 
the ratio w(t) = -jgA; the plot of the function w(t) is 
presented in the left panel of Fig. 12. Clearly, this func- 
tion is constrained: < w(r) < 1; maximal value of this 
ratio is about 1/3 in the first epoch of the Universe's evo- 
lution, when the DM can be considered as an effectively 
ultrarelativistic substrate. In the early Universe the DM 
state evolves quasiperiodically, i.e., the eras with effec- 
tive cooling were changed by the eras of effective heat- 
ing. Starting from t ~ 4 the function w(t) tends to zero 
linearly in t= log , i.e., the DM behaves effectively 
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W(t)=P/E 




FIG. 12. The left panel contains the plot of the ratio 
w(t)= , which presents the effective equation of state for 
the DM component of the dark fluid. Maximal value of this 
ratio is about 1/3 in the first epoch of the Universe's evo- 
lution, when the DM can be considered as an effectively ul- 
trarelativistic substrate. In the early Universe the DM state 
evolves quasiperiodically; starting from r ~ 4 the function 
w(t) tends to zero linearly in r, i.e., the DM behaves ef- 
fectively as a cold gas (dust). The right panel displays the 
cross-points of the p and E plots. This example illustrates 
the so-called coincidence problem: why the energy densities 
of the DE and DM components of the dark fluid are of the 
same order today (72% and 23%, respectively). 



as a cold gas (dust). 



C. Coincidence problem 

Although the DE and DM components of the dark fluid 
evolve at different rates throughout the history of the 
Universe, their magnitudes are of the same order today 
(72% and 23%, respectively) . This is known as the coin- 
cidence problem (3q|-(39|. The Archimedean-type model 
could help us to make a step toward solving this problem. 
In the right panel of Fig. 12 we placed the example of the 
plots describing the DE and DM evolution with time r. 
The guiding parameter p$ defines the asymptotical ratio 
between p and E; this ratio can in principle be chosen 
so that -§ 



|| . Other guiding parameters define how 
many cross-points of the p and E plots can exist. For in- 
stance, in Fig. 12 there are examples with one and three 
cross-points. Using the guiding parameters we can re- 
move one of the cross-points away from the point r ~ 1 , 
but we hope to discuss in detail the fitting problem in a 
special work. 

D. On the partition of the Universe's history into 
epochs and eras and its multi-inflationary behavior 

There are a few versions of dividing the history of 
the Universe into self-sufficient parts that are interest- 



ing from a physical point of view. For instance, one can 
consider the Hubble function H(t), find its zeros and ex- 
trema and separate the admissible time interval in line 
with them (see, e.g., [40]. We attached such a partition to 
the function —q(t), i.e., based the classification of mod- 
els on the zeros and extrema of the acceleration parame- 
ter. This seems to be motivated, since we are interested 
in picking out the epochs of accelerated and decelerated 
expansion of the Universe, thus the transition points ap- 
peared as points, in which the function — q(t) changes the 
sign. We have shown above that there are a lot of models 
in which the Universe's history is multistage. In partic- 
ular, the Universe's evolution can be quasiperiodic, and 
one can use the term multi-inflationary evolution, the 
first inflation being the sharpest, others being more and 
more smoothed. Following this line, we divide the epochs 
into eras by using the maximums and minimums of this 
function: this is equivalent to the method of statcfindcrs 
proposed in (34J. We hope to consider such fine details 
of the Universe's history partition in a special work. 



E. Conclusions 

The model of Archimedean-type coupling between 
dark energy and dark matter makes it possible to ex- 
plain the principal cornerstones of the Universe's evolu- 
tion: the early-time inflationary expansion, the late-time 
accelerated expansion, and the coincidence phenomenon. 

The model of Archimedean-type coupling between 
dark energy and dark matter possesses a wide set of 
guiding parameters suitable for fitting of the model pre- 
dictions to the observational data; we hope to devote a 
special work to this important question. 

The model of Archimedean-type coupling offers a nat- 
ural approach to the partition of the Universe's history 
into epochs and eras using the number of transition 
points and number of extrema of the function —q(t), the 
acceleration parameter of the Universe's evolution. The 
model allows us to speak about multistage evolution and 
about the multi-inflationary behavior of the Universe. 
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